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Abstract 

The Kubo formula is used to calculate the spin-Hall conductivity GsH in a spin-polarized two- 
dimensional electron system with Rashba- type spin-orbit interaction. As in the case of the unpo- 
larized electron system, GsH is entirely determined by states at the Fermi level, a property that 
maintains in the presence of magnetic impurities. In the clean limit, the spin-Hall conductivity 
decreases monotonically with the Zeeman splitting, a result of the ordering effect on the electron 
spins produced by the magnetic field. In the presence of magnetic impurities, the spin-dependent 
scattering determines a finite renormalization of the static part of the fully dressed vertex correc- 
tion of the velocity operator that leads to an enhancement of GsH-, an opposite behavior to that 
registered in the presence of spin-independent disorder. The variation of GsH with the strength of 
the Rashba coupling and the Zeeman splitting is studied. 

PACS numbers: 72.10.-d, 72.20.-1, 72.90.-fy 
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INTRODUCTION 



Equivalent to a local, momentum-dependent effective magnetic-field, the spin-orbit inter- 
action (SOI) in two-dimensional (2D) electron systems introduces a spin-dependent, chiral 
motion of the electrons that is sensitive to the application of an electric field. This property, 
that opens up the possibility of manipulating the electron spins exclusively by electrical 
means, is at the root of the tremendous amount of interest in understanding the electron 
dynamics in the presence of SOI, given the potential applications to spintronics. 

One such example is the intrinsic spin-Hall effect, when a pure spin current flows in a 
transverse direction under the action of an electric field The spin current is polar- 

ized along the third perpendicular direction. The magnitude of the spin-current response, 
described by the spin-Hall conductivity asH, reaches, in a clean system, a universal value 
e/Svr, independent of any sample parameters. 

The behavior of the intrinsic spin-Hall effect in the presence of non-magnetic impurities 
has been a subject of intense investigation. While analytic calculations led to a cancelation of 
the effect even in the presence of infinitesimal impurity concentration m , numerical studies, 
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done in finite size samples [2|, |5|, 16| , indicated that the spin-Hall effect persists in mesoscopic 
samples, up to a certain disorder strength. It has been shown that within the bulk, the spin- 
Hall conductivity is decaying exponentially along a distance of the order of magnitude of the 
spin precession length 7|]. More recent reports indicate that the discontinuous variation of 
asH in the infinite 2D system, from a finite value in the clean system to zero in the presence 
of the infinitesimal disorder, can be explained by introducing an additional dephasing effect 
associated with the inelastic electron lifetime |8|. This result suggests that the spin-Hall 
conductivity is enhanced by interactions that introduce additional scattering of the the 
electron spins and maintains a finite value even when the clean-disordered transition is 
performed. Naturally, one wonders if the opposite effect might be true. Are interactions 
leading to an ordering of the spins, such as the Zeeman coupling to an external magnetic 
field, acting as decreasing factors on cr^//? 

Inspired by these ideas, we proceed to a calculation of the spin Hall conductivity in a 
2D system with Rashba spin-orbit coupling, spin-polarized by a static magnetic field, per- 
pendicular on the sample. The alignment of the electron spins along the direction of the 
magnetic field counteracts the spatial disordering induced by the spin-orbit coupling, leading 



2 



in consequence to a diminished contribution to the spin current. Further, we consider mag- 
netic impurities and study the competing effects of the Zeeman sphtting and spin-dependent 
impurity scattering. The latter affects the magnitude of through the renormahzation 
effects it induces on the vertex corrections of the current operator. 

The simple model we discuss below, that of a non-interacting 2D spin-polarized electron 
gas with SOI and magnetic impurities, allows the simultaneous investigation of the intrinsic 
anomalous Hall effect, which would occur only when a finite magnetization is present, and 
of the spin-Hall effect in the presence of a distribution of magnetic scatterers, previously 
analyzed within a paramagnetic system ^. Our calculation is based on the Kubo formula, 
where we take into account the scattering of the electrons on the magnetic impurities. The 
algorithm discussed here, generalizes to spin transport the traditional treatment of the off- 
diagonal anomalous Hall conductivity of Ref. 1101, a method that has been also used with 
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great success to investigate the anomalous Hall effect in graphenes |12| . Within this 
framework we start by obtaining the impurity-averaged single-electron Green's functions 
and the renormalized vertex correction of the velocity operator. Then, we apply the Kubo 
formula to estimate the spin-Hall conductivity. First, in the case of a clean system, we use the 
exact eigenvalues-eigenstates of the Hamiltonian and obtain an analytic result for cTsh which 
show its dependence on the Zeeman splitting. Later, we use the impurity averaged Green's 
functions and the vertex-corrected current operator to estimate the spin-Hall conductivity in 
the presence of magnetic impurities. Analytical expressions for asH are derived as functions 
of the Zeeman splitting and the magnetic impurity scattering. 



THEORETICAL FRAMEWORK 



Model Hamiltonian 



We consider a non-interacting two-dimensional (2D) electron gas with Rashba-type spin- 
orbit coupling (proportional to the linear momentum) in the presence of a magnetic field. 
The system is assumed to contain magnetic impurities. The magnetic field B, is oriented 
along the z direction and is perpendicular on the layer. The resulting Zeeman splitting 
Ez = 275, proportional to the gyromagnetic factor 7 is considered a parameter of the 
problem. The noninteracting, single-particle Hamiltonian, written for an electron of wave- 
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vector k = {fcx; ^y} and kinetic energy = h'^k'^/2m in respect with the Fermi surface /i, 
(m is approximated by the bare mass) is 

Ho = ek + a{ky(T^ - k^dy) - Eza^ , (1) 

where a designates the spin-orbit couphng constant, while (Tj, (i = x, y, z) are the Pauh ma- 
trices. In the two-dimensional spin space, an elementary diagonalization procedure generates 
the two eigenvalues 



Sk,± = t,T + El (2) 

and the associated eigenstates of the Hamiltonian: 

I cosfe^<^/2 , , 
with cos 6* = Ezjl^k-, sin 6* = ak/^k and 





( sin |e^<^/2 \ 






I cos f 1 



Afc = y^«2^2 + El (4) 

the effective Rashba gap. 

When the magnetic impurities are present, an additional coupling Hamiltonian has to be 
included in Eq. ([1]): 

Himp = J s ■ S (5) 

with s and S denoting the electron and the impurity spin, respectively. The electron spin 
is treated like a quantum mechanical observable, described in terms of the spin-dependent 
creation and destruction operators at site i, Ci = Cii) and c| = (cj^, c]^) by s = |cj cr Cj. 
The impurity spin is considered to be a classical variable, whose direction n, in spherical 
coordinates is specified by the angles 9 and 0: S = Sn = S {sin 9 cos (p, sin 9 sin cj), cos 9). In 
the spinor representation, the coupling can be described by the matrix 

(cos6' sin6'e~*'^ \ 
\=(Tn. (6) 
sin — cos 9 I 

With u = hJ S/2, a rescaled exchange coupling, the interaction Hamiltonian is then written 



as 



Himp = M ^ c| (cr ■ n) d (7) 
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Throughout this analysis, the impurity scattering problem is treated perturbatively, as we 
neglect the regime where the Kondo effect may be important. In our approximation, the 
lifetime of the quasiparticles at the Fermi level is evaluated for each band, as the imaginary 
part of the self-energy in the second order perturbation theory. At the same time, the shift 
of the chemical potential, due to the real part of the self-energy, is not considered. 



Green's function, self-energy and current vertex correction 

The free electron Green's function is obtained from the single particle Hamiltonian, 
Eq. ([1]) as a 2 X 2 matrix in the spin space: 

u -ek + fi + a{kyO.^ - k^cXy) - Ezcr^ 



(8) 



[lu — £k,+ + fi + i6 sgn (u)] [cu — £k- + fi + i6 sgn{uj)] 

with (5 > an infinitesimally small quantity. In the presence of the impurities, G^{uj) is 
modified to include the effects of the elastic scattering. The relaxation time is given by the 
imaginary part of the self-energy, which, in the lowest order (see Fig. [T]), is obtained from: 



(9) 



As an explicit function of the effective Rashba gap and the Zeeman splitting, the self-energy 



FIG. 1: Second order self-energy contribution due to magnetic impurities. The solid line represents 
the Green's function while the crosses describe the interaction with a single magnetic impurity. 



is given by: 



X 



(27r)2 J A-K 
1 



cr ■ n 



uj — + + ^ + i5sgn{uj) — — A^. -|- /i -|- i5sgn{uj) 



cr ■ n 



(10) 
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cr • n 



uj — tk + + ^^ + iSsgn{uj) cu — ek — + 1^ + i5sgn(uj) 
Since the real part of Eq. (fTOj) just renormalizes the Fermi energy, we focus only on its 
imaginary part, the one that determines the quasiparticle lifetime at the Fermi level. In 



contrast to the case of non-magnetic impurities, now the scattering rates depend on the 
chirahty of the band: 

4:71 



-"^^"(^^ 2 J {2ny 
(T ■ n) [6{uj - ek + Ak + fi) + 6{uj - ek - Ak + fi)] {cr ■ n) 



(11) 



+ cr ■ n 



A. 



[6{uj -ek + Ak + /i) - S{uj -ek- Ak + /i)] (cr • n) 



The momentum space integral is processed by changing to an integral over energy, 
J d?k/(2'KY J deNo{e) where A'o(e) = m/2iT the density of states at the Fermi surface. 
After performing the integrals over the solid angles, and some standard manipulations, we 
finally obtain: 



We recognize that 



m 



1 - 



+ 0{fi - Ez) 



1 - 



(12) 



1 



Afc 



F+ 



(13) 



are the densities of states in the chiral bands, allowing us to define the symmetric and 
antisymmetric scattering rates: 



N+ 



where Eq. (jl]) was employed. Thus, 

QmHiiuj) 



1 (\ 1 



„ , -^-Oz] sgnl^wj 

I \ T T 



(14) 



(15) 



We introduce the band-dependent impurity scattering times, r^^ = ± and define 
r = l/2r, r = l/2f, r± = 1/2t± . with these notations, the impurity-averaged Green's 
function G'fc(e) becomes 

u + ir sgn(u;) - + /i + a{kya^ - k^Oy) - [Ez + if sg\i{uj)]az 



Gk{uj) 



(16) 



[ijj - 8k,+ + + «r+ sgn(cj)] [uj - 8k- + /i + iT_ sgn(co')] 

Eq. flTB]) generates the retarded (R) and advanced (A) Green's functions, G^j^'^\uj), in the 
second-order perturbation theory. 

The next ingredient needed for computing the spin-Hall conductivity is the current vertex, 
involved in the calculation of the polarization bubble when multiple scattering events on the 



FIG. 2: Diagrammatic representation of the current vertex. The dots represents the bare current, 
the upper/lower soHd Hues are the retarded/advanced Green's function and the crosses represent 
the interaction with the magnetic impurities. 

magnetic impurities are considered. We determine the vertex-corrected current as the ladder 
series expressed in Fig. |2l whose equivalent analytical equation is: 

47 J (2^ ^'^ ■ ^^^'^ "^^^^^^ + ^){^- n) (17) 
In the static limit, when e — > and ^ 0, we write 

7. = ^ - + J ^{^- n)Gf (0) 7. G^(0)(^ ■ n) (18) 

where the advanced (A) and retarded (R) electron Green's functions are obtained from the 
static limit of Eq. ( |T6l) . A solution to Eq. ( |T8l) can be obtained in the form of an expansion: 
1x = kx/m — 'JxCfj,. First, a simple analysis shows that two components of the static part of 
the dressed vertex vanish: 7^ = 7;^ = 0. The remaining, non-zero component of the vertex 
function is 7^ expressed as: 

y = (] ^ f tkjtk - At) \ / n, C (Pk (efc - ^Jif - El 

" " V 3 J (27r)2 A+A. )[ ^ 3 J {2ny A+A. 

where 

/i± = (/i-^,±)Vr|. (20) 

In the case of non-magnetic impurities the vertex coefficient 7^ cancels when the Fermi level 
is in the upper band (both bands are occupied), leading to the disappearance of the spin-Hall 
effect in the thermodynamic limit when any amount of disorder is present. This is not the 
case, however, when magnetic impurity scattering occurs, since now the static part of the 
dressed vertex is larger than the bare Rashba coupling for any ratio fi/Ez > — 1. In our 
model the Rashba coupling is the static part of the bare vertex. 

As can be seen in Fig. |3l in the extreme case when the chemical potential satisfies 
= Ez + 5e (with 5e some positive infinitesimal energy), so just one band is partially 




1 2 3 4 5 6 

a (10 "^^eVm) 

FIG. 3: jx as function of the Rashba coupling strength for different ratios of chemical potential 
versus the magnetic energy fi/Ez- Here Ez = 1.1 x 10~'^eV and the effective electron mass is the 
bare one me- 

occupied, the vertex is practically not renormalized and takes the bare value a. The renor- 
malization is more pronounced as the band is gradually filled. The largest renormalization 
is obtained when both bands are occupied /i > Ez + Se ■ Increasing the ratio fi/Ez above 1 
does not lead to a larger renormalization. Therefore, for any filling factor, in the thermody- 
namic limit, the spin-Hall conductivity is finite when magnetic impurities are present in the 
system, irrespective of how strong/weak the interaction potential is. For the experimentally 
accessible values of the Rashba coupling strength of 5 — 6 x [13] we present the 

behavior of the vertex function in Fig. [3l An external magnetic field of 5 = IT is considered 
to be applied, so that the Zeeman energy is approximately Ez = 1.1 x 10~'^eV. 

SPIN-HALL CONDUCTIVITY 

In this section we present analytical results for the spin-Hall conductivity. First we 
compute for a non-interacting system. Because in the clean limit the eigenvectors and 
eigenenergies are exactly known, the Kubo formula can be employed, and an analytical result 
is obtained in agreement with previous analytical work. When magnetic impurity effects 
are investigated, Kubo formalism in term of exact eigenstates/eigenenergies is no longer 
suited and the causal Green's function method has to be considered. For that, the impurity 
averaged Green's functions and the vertex correction derived in the previous section are 
needed. 
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Spin Hall conductivity for the non-interacting system. Exact result 



The Kubo formula that determines the spin-Hall conductivity in the clean system is 
written in the chiral basis of states, Eq. 



(i^k „ < kn'\il\kn X kn\vJkn' > 
Ukw - /fen) 3m ^ 



(2vr) 



kn' 



(21) 



where n is a band index, in our case n, n' = ±. The electron velocity along the y direction 
is Vy = ky/m + aax and the i-polarized current propagating in the x direction is H = 
h/4{vx, (Jz} = hkx/2maz. Upon the insertion of their matrix elements, evaluated in the 
chiral basis, in Eq. fl2Tl) we obtain: 

k^dk 



a. 



xy 



IGvrm 



-{fk+ — fk~) 



(22) 



with fk± the Fermi-Dirac distributions corresponding to the two bands. In the absence of 
the magnetic field, when Ez — > 0, we recover the well known result 1| a^y = e/Svr for a 
clean two dimensional electronic system. For a finite Zeeman splitting, an analytical result 
can be derived in terms of the Fermi energies of the chiral bands: 



xy 



in ma' 



El 



ma tF+ 



E\ + ma^ep_ 



(23) 



+ 2ma2eF+ \/ E\ + Ima^ep- 

a result that shows that even the simple presence of an external magnetic field leads to a 
non-universal value for the spin-Hall conductivity. 



Spin-Hall conductivity using the causal Green's function. The role of disorder 

In this section we derive an analytical expression for the spin-Hall conductivity when 
both magnetic impurities and external magnetic field are considered. The Kubo formula for 
the spin-Hall conductivity written for the impurity averaged Green's function gives: 

-^.M = ^^'1^1 ^Afyfi^)iG^i^)-G'{e))lxG\e-u) 

-fyG^hJie - c.)(G«(6 -co)- G\e - co))) (24) 

where the Green's functions include the scattering lifetimes and < . . . > represents the im- 
purity configuration average. There are two types of contributions to the integral, Eq. (12^ : 
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a = 0.5xlO"''eVm 

■■■■ a = 1.0x10"'' eVm 

-- a = 3.0xlO"''eVm 

-- a = 5.0xlO"''eVm 



0.5 1 

Magnetic Field ( T ) 



1.5 



FIG. 4: Magnetic field dependence of the spin-Hall conductivity for different SOI couplings when 
no vertex correction is considered (7I a). For the case of magnetic impurities this is a good 
approximation. Here fJ-fEz = 1.5 

from states below the Fermi level and from states close to the Fermi level [l4|. The con- 
tribution from states well below the Fermi level can be neglected in the limit app <^ cf, 
or eFT± ^ 1, because it contains only combinations of the form G^G^ and G'^G^ [15]. At 
the same time, magnetic impurities, have practically no effects on these states due to small 
scattering rates, as compared to their energies. In stark contrast, states close to the Fermi 
energy are strongly effected once their energy becomes comparable to the scattering time 
l/r±. 

The vertex correction becomes important when processes across the Fermi surface are 
considered, so that one G^^^ and one G^"^^ enters the vertex equation ([T8|) . By taking the 




a (10 eV m) 



FIG. 5: Spin-Hall conductivity as function of the Rashba spin-orbit interaction strength for different 
ratios n/Ez- Here an external field of B = IT is considered. 
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zero frequency limit in Eq. (12^ the spin-Hall conductivity becomes: 

'''' 

with A± given by Eq. (!20l) . Simple algebraic manipulation gives an exact result for the 
spin-Hall conductivity in terms of the density of states for the chiral bands. 

= ^ 4™,jr^j-l^ (iV„(f0| + A'oCf^) (26) 

where we have introduced the quantities: 81^2 = + ± {m^a^ + 2jjLma^ + and 

/i - gi,2 - (^1 + 2ma^£,^2) /i - ^1,2 + (^| + 2ma^gi,2) . „^ 

^ 2r+(^i-f2) ^ 2r_(^i-^2) ^ ' 

In the absence of magnetization {Ez — > 0), when no vertex corrections are considered 
(7^ — i> a), and in the weak disorder limit, Eq. (!25|) generates the well-known universal 
expression of the spin Hall conductivity, a^y = e/Sir. One important observation is that, 
in contrast to the case of unpolarized disorder, when the static component of the dressed 
velocity is renormalized to zero by the vertex correction, here the vertex corrections lead 
to an enhancement of 20 — 50% of the bare static velocity (see Fig. [2])- This observations 
shows that, in the case of magnetic impurities, even when vertex correction for the velocity 
are neglected, a good enough approximation for the spin-Hall conductivity is obtained. In 
this limit (7^ — >■ a) we present in Fig. H] and Fig. [5] typical behaviors for the spin-Hall 
conductivity as function of the external field as well as function of the spin-orbit strength 
(similar curves are obtained also when the vertex is considered). Typically, for a given a, 
the magnetic field reduces the strength of the spin-Hall effect. This can be, in principle 
understood, by considering the different polarization effects of the magnetic field, that stati- 
cally orientates the electron spins from both chiral bands along its direction, and the Rashba 
interaction that induces an in-plane, dynamic polarization. The larger the ratio ma^ / Ez is, 
the stronger the spin-Hall effect is, and in limit of zero magnetic field the universal expression 
for the spin-Hall conductivity is reobtained. 



CONCLUSIONS 



The present work addresses an important topic in the field of spin-Hall effect, that is the 
effect of magnetic impurities and the role of a Zeeman term on the spin-Hall conductivity. 
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We have obtained simple but robust analytical results for the spin-Hall conductivity, which 
in some particular limits converge to the previous known results ^, ^• 

First we find that the spin-Hall conductivity is no longer universal in the presence of a 
magnetic field even in the clean limit. The most important observation is related to the 
behavior of the system when magnetic impurities are present. In this case vertex correction 
leads to an enhance of the spin-Hall effect, contrary to the case of non-magnetic impurities 
where the static part of the fully dressed vertex identically vanishes in the weak scattering 
limit. This allows us to conclude that the bare vertex is a good approximation when com- 
puting the spin-Hall conductivity, and that the bare bubble diagram is good enough when 
computing the spin-Hall conductivity in the presence of magnetic impurities. 
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